Deformation of the N=Z nucleus 76 Sr using /3-decay studies 
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A novel method of deducing the deformation of the N=Z nucleus 76 Sr is presented. It is based on 
the comparison of the experimental Gamow- Teller strength distribution B(GT) from its /3-decay with 
the results of QRPA calculations. This method confirms previous indications of the strong prolate 
deformation of this nucleus in a totally independent way. The measurement has been carried out 
with a large Total Absorption gamma Spectrometer, "Lucrecia" , newly installed at CERN-ISOLDE. 
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The shape of the atomic nucleus is conceptually one 
of the simplest of its macroscopic properties to visualise. 
However, it turns out to be one of the more difficult prop- 
erties to measure. In general terms we now have a picture 
of how the nuclear shape varies across the Segre Chart. 
Nuclei near to the closed shells are spherical. In con- 
trast nuclei with the valence nucleons in between two 
shells have deformed shapes with axial symmetry and 
the extent of the quadrupole deformation is quite well 
described as being proportional to the product N p N„ of 
the numbers of pairs of valence protons (N p ) and neu- 
trons (N„) Q. This picture is underpinned by both the 
Shell and Mean Field models of nuclear structure. Exper- 
iment and theory concur that, as the N p N„ parameter- 
isation would suggest, nuclei rapidly deform as we add 
only a small number of valence nucleons to the magic 
numbers. Thus nuclei in the middle of the £7/2 shell turn 
out to be deformed even although the numbers of valence 
nucleons are relatively small. 

Experimentally this picture is supported by a mass 
of independent observations: the strongly enhanced 
quadrupole transition rates between low-lying states, the 
strongly developed rotational bands built on low-lying 
states, and measurements of ground state quadrupole 
moments. Where we have evidence of the shapes of 
ground and excited states in the same nucleus they are, 
in general but not always, the same. It turns out that in 
some cases nuclear states with different shapes co-exist 
in the same nucleus 0. 

The nuclei with NwZ and A«70-80 are of particular 
interest in this context. Such nuclei enjoy a particular 



symmetry since the neutrons and protons are filling the 
same orbits. This, and the low single-particle level den- 
sity, lead to rapid changes in deformation with the ad- 
dition or subtraction of only a few nucleons. In terms 
of Mean Field models these rapid changes arise because 
of the proximity in energy of large energy gaps for pro- 
tons and neutrons at Z,N=34 and 36 on the oblate and 
Z,N=38 on the prolate side of the Nilsson diagram. As 
a result Mean Field calculations predict the existence of 
several energy minima with quite different shapes in some 
of these nuclei 0, • Evidence of this co-existence has 
been found for instance in Se and Kr nuclei laJfj, and 
it is also predicted for the lightest Sr isotopes 0. Thus 
it is of considerable interest to map out the deformation 
of the ground and excited states of nuclei in this region. 
This is easier said than done, however. There are a num- 
ber of methods to measure the deformation of the ground 
state in unstable nuclei based on the interaction of the 
electric quadrupole moment of the nucleus with an exter- 
nal electric field gradient || . These techniques are not 
applicable to nuclei with J=0 or 1/2, moreover they very 
seldom give the sign of the quadrupole moment and hence 
cannot distinguish between oblate and prolate shapes. 

Here we present an alternative method to deduce 
whether the ground state shape of an unstable nucleus is 
oblate or prolate, and apply it to the N=Z=38 nucleus 
76 Sr, the most deformed nucleus in the region accord- 
ing to Mean Field calculations and previous in-beam 
experiments 01 • The method is based on an accurate 
measurement of the Gamow- Teller strength distribution, 
B(GT), as a function of excitation energy in the daughter 
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nucleus, and relies on the technique of Total Absorption 
gamma Spectroscopy (TAgS) which will be explained 
later. The theoretical idea was suggested b y H amamoto 
et al. ^3 an d pursued by Sarriguren et al. |l2| . Accord- 
ing to them, one can study the deformation of the ground 
state of a particular nucleus by measuring the B(GT) dis- 
tribution of its /3-decay. In these references the authors 
calculate the B(GT) distributions for various nuclei in 
the region for the deformations minimising the ground 
state energy. In some cases, the results differ markedly 
with the shape of the ground state of the parent, espe- 
cially for the light Kr and Sr isotopes. 

A precise determination of the B(GT) distribution is 
required for such studies and this is far from trivial. Tra- 
ditional high resolution techniques, based on the use of 
high purity Germanium (HPGe) detectors to measure the 
7-rays emitted after the /3-decay, often fail to detect sig- 
nificant but very fragmented strength at high excitation 
energy in the daughter nucleus. This is mainly due to 
three factors: the low photo-peak efficiency of HPGe de- 
tectors for high energy 7-rays, the high fragmentation of 
the B(GT) at high excitation energy, and the fragmen- 
tation of the gamma de-excitation of the levels in the 
daughter through many different gamma cascades. To- 
gether they cause the so-called Pandemonium effect |l3l | : 
many weak cascades de-exciting levels at high energy can 
remain undetected leading to large systematic errors in 
the determination of the B(GT). This is the reason why, 
even although Refs. 0] and 0] give the first indication 
of the prolate character of the 76 Sr ground state, one 
must determine the B(GT) distribution more accurately 
over the whole Qec window to provide conclusive proof. 

The alternative, the TAgS technique, avoids these sys- 
tematic uncertainties. The basis of this method is the de- 
tection of the entire energy of the gamma cascades rather 
than individual 7-rays. For this purpose one needs a high 
efficiency detector with acceptable resolution for gammas 
such as the inorganic scintillators Nal(Tl) or BaF2. Fur- 
thermore, this detector must have a geometry as close 
as possible to 4ir to absorb the complete cascade energy. 
If this is achieved one can measure directly the /3 in- 
tensity 1/3 (E) as a function of the excitation energy in 
the daughter nucleus, and from this one can extract the 
B(GT) distribution. 

In this paper we present the results of TAgS measure- 
ments on the decay of the N=Z nucleus 76 Sr. A com- 
parison of the results of this measurement with the cal- 
culations of Ref. 01 allows us to establish the prolate 
character of the ground state of 76 Sr without ambiguity. 

With the aim of measuring the /3-decay of nuclei far 
away from the stability line with the Total Absorption 
technique, a spectrometer called "Lucrecia" has been 
installed at the ISOLDE mass separator at CERN. It 
consists of a large Nal(Tl) crystal of cylindrical shape 
(L=0=38 cm) with a cylindrical hole (0=7.5 cm) at 
right angles to the symmetry axis. The purpose of the 



hole is twofold: on the one hand it allows the beam pipe 
(coming from the separator) to enter up to the centre of 
the crystal, thus allowing on-line activity of very short 
half-life (>5 ms) to be deposited here and measured. On 
the other hand it allows us to place ancillary detectors 
inside for the detection of the positrons (/3 + -decay), elec- 
trons (/3~ -decay) or X-rays (EC process) produced in the 
decay. Surrounding the whole setup there is a shielding 
box 19 cm thick made of four layers: polyethylene-lead- 
copper-aluminium (see Ref. for further details). 

In order to produce the nucleus of interest ( 76 Sr), a 
52 g/cm 2 Nb target was bombarded with a 1.4 GeV pro- 
ton beam. The intensity was chosen to produce a count- 
ing rate of «3.5 kHz in the Nal crystal. In order to sepa- 
rate Sr selectively, a fluorination technique was used [ljj • 
The radioactive beam was steered to the detector setup 
and implanted in an aluminised mylar tape which was 
moved every 15 seconds to transport the source to the 
middle of the crystal and to avoid the buildup of the 
daughter activity (T 1/2 ( 76 Sr)=8.9 s, T 1/2 ( 76 Rb)=36.8 s). 
During this 15 s cycle the decay of the implanted ra- 
dioactive source was measured. The 7-rays following the 
decay (either by (3 + or by EC) were measured by the 
Nal(Tl) crystal and analysed without any condition on 
the ancillary detectors. However these detectors were 
very useful for the on-line control of the measurement. 

In ideal conditions, if the TAgS had 100% peak ef- 
ficiency over the whole energy range, the experimental 
spectrum measured in the Nal(Tl) cylinder would be the 
j3 intensity distribution Ig(E) convoluted with the energy 
resolution of the crystal and the response of the detector 
to the positron when applicable. In reality the detector 
does not have 100% peak efficiency because of the dead 
material inside the spectrometer (the ancillary detectors) 
and the transverse hole. Consequently the spectrum is 
modified by the response function of the detector. In 
other words, the relationship between the quantity of in- 
terest, I/3 (E), and the experimental data d(i) is: 

z — J \j = energy bin J 

In order to obtain Ia(E) from our data we should solve 
Eq. QJ. This is not a trivial task because the response 
matrix R(i,j) can not be inverted due to the fact that 
it is quasi-singular in the sense that two neighbouring 
columns are very similar. However, there is a set of al- 
gorithms that has been developed to solve this kind of 
"111 Posed" problem. In Ref. [18j there is a systematic 
study of three of these methods applied to the specific 
problem of the TAgS data. Here we have used the Expec- 
tation Maximisation algorithm [I<| to obtain the L3 (E) 
by unfolding the experimental data. To calculate the re- 
sponse matrix R(i,j), which is needed by the algorithm, 
we have used the levels and branching ratios given in 
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Ref. |15|. and the GEANT4 simulation code. The analy- 
sis has been performed taking into account both the EC 
and /3 + components of the decay. A more detailed ex- 
planation on the procedure to calculate R(i,j) and to 
analyse the data will be given in a forthcoming article. 

The best check one can perform to validate the result 
of the analysis is to recalculate the experimental spec- 
trum by multiplying the response function of the detector 
(R(i,j) in Eq. JJJ) by the resulting beta intensity lp. If 
the analysis is properly done, this recalculated spectrum 
should be very similar to the real experimental spectrum. 
The upper part of Fig. ^ shows the experimental spec- 
trum (shade without line) overlaid with the recalculated 
one (dashed line). The agreement between the two spec- 
tra is very good. 

The 1/3 (E) is the experimental result of this work, how- 
ever, the physical information is carried by the reduced 
transition probability B(GT), which can be extracted 
from the I,g(E) using the expression: 
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B{GT) = -—IA- x614 7 (^y (2) 
J{Qec ~ E)T 1/2 \g A J 

where the B(GT) is averaged inside the 40 keV energy 
bin, and /(Qec — E) is the Fermi integral which carries 
the information on both the phase space available in the 
final state and the Coulomb interaction. For the calcu- 
lation of the B(GT) we have used the Qec value from 
Ref. [U UK, the T 1/2 from ^ and the tabulated Fermi 
integral from (2^ . 

In the lower panel of Fig.nthe resulting B(GT) distri- 
bution is presented. The analysis gives a total B(GT) of 
3.8(6)g^/47r up to 5.6 MeV of which 57% is located in the 
resonance between 4 and 5 MeV. This resonance is weakly 
visible in the almost structureless TAgS spectrum of the 
upper panel. Its large B(GT) value is a consequence of 
the strong dependence of the Fermi integral with the en- 
ergy. At lower energy, the marked B(GT) to levels at 
0.5 MeV, 1.0 MeV and 2.1 MeV were already observed 
in 0|, although the B(GT) values are in disagreement 
with our results due to the already mentioned Pandemo- 
nium effect. It should be noted that the /3-delayed proton 
emission (S p =3.5 MeV) in this decay has been observed 
at excitation energies from 4.8 to 5.8 MeV 0,^3- How- 
ever, the contribution of this component is very small, of 
the order of 2% in B(GT) compared to the decay through 
the /3-delayed 7-rays studied here. 

One way to compare the results with the theory is to 
accumulate in each energy bin the sum of the B(GT) 
measured up to that particular energy. Fig. [2] shows this 
plot in which the experimental result is compared with 
the theoretical calculations of Ref. 0] for both pure pro- 
late and oblate shapes for the ground state of 76 Sr. The 
generally accepted quenching factor of 0.6 has been ap- 
plied to the calculations. 



Figure 1: Upper panel: Experimental total absorption spec- 
trum of the /3-decay of 76 Sr overlaid with the recalculated 
spectrum after the analysis (see text). Lower panel: B(GT) 
distribution derived from the experimental data shown above. 
The shade represents the experimental uncertainty. 

In Ref. ^3 the authors first construct the quasiparti- 
cle basis self-consistently from a deformed Hartree-Fock 
(HF) calculation with density-dependent Skyrme forces 
and pairing correlations in the BCS framework. The 
minima in the total HF energy vs deformation param- 
eter plot give the possible deformations of the ground 
state. For the case of 76 Sr two minima are found, 
one prolate with /?2=0.41 and the other oblate with 
/32=-0.13. Finally, the QRPA equations are solved with 
a separable residual interaction derived from the same 
Skyrme force used in the HF calculation. For the B(GT) 
calculation the same deformation is assumed for the 
ground state of the parent and for the levels populated 
in the daughter nucleus. Fig. |2] shows the results using 
the residual interaction SK3. The agreement of the 
experimental results of this work (s quar es in Fig. |2J 
with the prolate shape calculation of [12( is very good 
over the energy range 0-5.6 MeV. In contrast, there is no 
similarity between the results of the oblate calculation 
and the experimental points. 

This agrees with the strong deformation (/?2~0.4) of 
76 Sr already extracted from the dynamical properties 
observed in in-beam studies |10|. It also gives the first 
definitive experimental evidence of the prolate character 
of the ground state deformation, confirming the result 
indicated in 0] and [lfij . 

We should point out here that these results prove 
the validity of the method of deducing the sign of the 
electric quadrupole moment of ground states or /3-decay 
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and with previous 
13. An important 
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Figure 2: Accumulated B(GT) as a function of the excita- 
tion energy in the daughter nucleus. The experimental re- 
sults from this work (squares) are compared with the theoret- 
ical calculations f 2] assuming prolate (solid line) and oblate 
(dashed line) shapes for the 76 Sr ground state. The shade 
indicates the experimental uncertainty. 



isomers from the study of their decay. This opens new 
opportunities in the study of nuclei far from the stability 
line where very often the first information comes from 
/3-decay (half- life, J"...). On the other hand the theo- 
retical approach used in the present study has been so 
far restricted to nuclei in this region. The present work 
should encourage further theoretical studies in other 
regions of well deformed nuclei. It is worth noting here 
that QRPA calculations have been successfully applied 
to describe other properties in nuclei with open shells 
where the nucleon pairing correlations are imp ortant 
such as B(E2) values [2^] or giant resonances |24| . 

Finally, as a part of this series of experiments, we 
have recently published the results on 74 Kr decay [lij . 
In that paper a clear indication of shape mixing was 
deduced, a conclusion which is further corroborated by 
the present study which is free of shape admixtures. 

Summarising, in this work we present the results 
of an experiment devoted to measuring the B(GT) 
distribution in the decay of the N=Z isotope 76 Sr. When 
we compare our experimental results with the theoretical 
calculations of Ref. 01 we conclude that the ground 



consequence of the present work is the validation of the 
method of deducing the deformation, including the sign 
of the quadrupole moment, from the comparison of the 
/3-decay TAgS results and the calculated B(GT) since 
the 76 Sr ground state is a very clean case, free of shape 
admixtures. 
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